Random matrix theory for closed quantum dots with weak spin-orbit coupling 
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To lowest order in the coupling strength, the spin-orbit coupling in quantum dots results in a 
spin-dependent Aharonov-Bohm flux. This flux decouples the spin-up and -down random matrix 
theory ensembles of the quantum dot. We employ this ensemble and find significant changes in the 
distribution of the Coulomb blockade peak height, in particular a decrease of the width of the distri- 
bution. The puzzling disagreement between standard random matrix theory and the experimental 
distributions by Patel et al. might possibly be attributed to these spin-orbit effects. 
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The spin-orbit coupling in a two-dimensional semi- 
conductor quantum well mainly contributes through the 
Rashba lj and Dresselhaus |2] terms, arising from the 
asymmetry of the confining potential and the lattice 
structure, respectively. It is much weaker than in three- 
dimensional semiconductors where it is induced mainly 
by impurities, which are absent in a high-mobility two- 
dimensional electron gas. The spin-orbit scattering is 
further suppressed if the two-dimensional system is con- 
fined to a quantum dot; estimates of the spin-flip rates 
were given m 

Ref.0 

This fact is of great importance for 
future applications of quantum dots as spintronics de- 
vices. However, it was shown that spin-orbit scattering 
has a significant effect in the presence of an in-plane mag- 
netic field 0- FJl ■ which explains Q recent experiments 



In this paper, we discuss another manifestation of the 
spin-orbit coupling in confined structures, which takes 
place even in the absence of appreciable spin-flip scatter- 
ing. Aleiner and Fal'ko recently showed |j| that a weak 
spin-orbit coupling creates a spin-dependent Aharonov- 
Bohm flux. While this flux does not flip spins, it can 
change the random matrix ensemble of the quantum dot. 
For broken time reversal symmetry, the spin-up and spin- 
down parts of the spectrum are completely uncorrelated 
and described by independent Gaussian unitary ensem- 
bles (GUE) 0. The possibility of such an ensemble was 
raised by Alhassid |q, while the relation to the spin- 
orbit coupling was already suggested by Lyanda-Geller 
and Mirlin |9(. In the present paper, we study the statis- 
tical distribution of the Coulomb blockade peak height in 
this ensemble, and find the distribution to be narrowed. 
This might explain the discrepancy between a recent ex- 
periment by Patel et al. [lfj and standard random matrix 
theory (RMT) dEl at low temperatures. 

In Ref. the free-electron Hamiltonian with Rashba 
and Dresselhaus spin-orbit terms was expanded to sec- 
ond order in the coordinates, under the assumption that 
L1.2/M.2 *C 1 (Li,2- lateral dimensions of the two- 
dimensional quantum dot; Ai^: characteristic length 
scale of the spin-orbit coupling which is proportional to 



the inverse spin-orbit coupling strength). One obtains 
H = ^(p-eA-a ± ^-a {l ) 2 +u(r). (1) 

Here, u(f) is the (disordered) confining potential; p = 
P — eA is the kinetic momentum with the canonical mo- 
mentum P and the vector potential 

A = B z [rxn z ]/2c; a ± = [rxn z ]/ (2\ 1 \ 2 ); (2) 
1 \rxn z ] ( x\<j\ x 2 <J2 



a 



6 A1A2 V ^1 



A 2 



(3) 



o~i denote the Pauli matrices and B is the magnetic field 
in the direction [001] perpendicular to the lateral quan- 
tum dot. The coordinates x\ and X2 are along the di- 
rections [110] and [110] and we neglected the Zeeman 
splitting term as we are interested in the behavior at rel- 
atively low magnetic fields. The term an is responsible 
for spin-flips but it is of higher order in the spin-orbit 
coupling strength than aj_. Thus, it will be neglected in 
the following as we assume the spin-orbit coupling to be 
weak such that a± dominates. The aj_ term has exactly 
the same form as the vector potential A except for its 
spin-dependence. As an electron collects an Aharonov- 
Bohm flux on a close path due to the vector potential 
A, it also collects a spin-dependent flux due to a±. This 
spin dependent flux translates to a spin-dependent ef- 
fective magnetic field, so that the electrons feel a to- 
tal magnetic field of strength B° s = B + ~ x~\j§ with 
(7 = ±h for up- and down-spin, respectively. An increase 
of the flux changes the matrix elements, and scrambles 
the eigenenergies and eigenvectors. In the absence of 
spin-orbit coupling, the flux is exactly the same for spin- 
up and spin-down electrons such that their eigenenergies 
and eigenvectors are degenerate. If the spin-orbit terms 
are present, but no external magnetic field is applied, the 
time-reversal symmetry is preserved, and the states are 
still Kramers degenerate (up-spin and down-spin see the 
same magnitude of magnetic field with opposite signs). 
However, when spin-orbit coupling and external magnetic 
field are present, electrons with different spin see different 
magnetic fields, and their eigenenergies and eigenvectors 
decorrelate. If the spin-dependent flux is large enough 



2 



O N=100 
X N=200 
□ N=400 
/ N=800 
- N=1600 




FIG. 1: Correlations of the level tunneling rates Cr and the 
rescaled spectral diffusion correlator Ce for the RMT model 
jH, as a function of x\/~N for different matrix sizes N. The 
data collapse for different N indicates the universality. 

spin-up and spin-down eigenenergies and eigenvectors are 
distributed according to two independent GUEs • 

Before we analyze this weak spin-orbit RMT ensem- 
ble, we study the decorrelation of the eigenenergies and 
eigenvectors due to a change in the magnetic field, to de- 
termine how much flux is needed in order to have two 
uncorrelated ensembles. The additional flux translates 
to a change of the random matrix described by |l2| 



H = 



Hi + xH 2 



(4) 



with RMT matrices Hi and H 2 in the unitary ensemble. 
As the perturbation x increases from zero the eigenen- 
ergies Ei(x) and eigenfunctions ^%{x) of H change. We 
analyze the decorrelations of the energies via the level dif- 
fusion correlator Ce = ((\/ (Ei(x) — Ei(0)) 2 )) / A, where 
((••■)) means averaging over different realizations and dif- 
ferent levels i. This correlator has been shown to have a 
universal form |13| . The decorrelation of the eigenfunc- 
tions is measured by Cr = ((| (tpi(x)\i(>i(Q)) | 2 )). It can 
be shown that this correlator also measures the corre- 
lations of the level tunneling rates and has a universal 
form as well The results are presented in Fig.^ an d 
show that the correlations in both quantities disappear at 
about the same value of xyfN « 1, where N is the size of 
the random matrix. Hence we conclude that the decorre- 
lation of the eigenvalues and the eigenfunctions (dot-lead 
coupling) occur together. Thus, the spin-orbit coupling 
leads to a crossover from two degenerate GUE spectra, 
to an ensemble of two uncorrelated GUE spectra. 

For the above RMT model, the crossover to two in- 
dependent GUE ensembles occurs at x\/N ~ 1. The 
corresponding flux difference needed to decorrelate the 
spectrum is given by the following relation |l2| 



xVN = xVgr 



g&eff 



(5) 



where <5$ c ff is the flux difference, <&o is the quantum unit 
of flux, gx denotes the Thouless conductance, and x is a 



non-universal sample-dependent constant of order unity. 
We thus realize that one needs about 1/ sfgr flux quanta 
to crossover to two uncorrelated GUE ensembles. 

Let us now estimate the strength of spin-orbit inter- 
action required to create this amount of flux difference. 
As mentioned above, the difference in effective flux be- 
tween the two spin sectors is S$ c ff/$o = LiL 2 /(XiX 2 ). 
The A's are connected to the Rashba and Dresselhaus 
spin orbit parameters 7 and r\ via I/A1A2 = 4(7 2 — n 2 ). 
Independent estimates of 7 and r\ are not available, 
but, in principle, can be obtained 0. One can get 
an approximate value via the better known parameter 
Q 2 = (hv F /E F ) 2 ( 1 2 +r 1 2 ): 
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Typical experimental values are k F L\^ 2 ~ 50, gx ~ 10 — 
100, and estimates for Qso are in the range 4 — 16 x 10 
0- Thus, we estimate ^fg~T~ 5$ e ff/$ = 0.1 - 6.4, i.e., 
the right hand side of Eq. (0) can be expected to be of 
order unity [l5| . Hence, the spin-orbit effect is strong 
enough to decorrelate (or to start to decorrelate) the 
spin-up and spin-down sector, while being weak enough 
not to yield significant spin-scattering. A strong spin- 
scattering, which can be generated by the application of 
an in-plane magnetic field, would mix the two spin species 
and result in a single GUE. 

We will now analyze this situation where the weak 
spin-orbit coupling results in two uncorrelated GUE en- 
sembles for spin-up and spin-down electrons and the re- 
sults do not depend on the spin-orbit strength. Under 
the assumption of a constant Coulomb interaction [l^ 
and applying the Master equation for sequential tunnel- 
ing through the quantum dot, the conductance of a closed 
quantum dot is given by [lq | (for a review see Ref. |g) 
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Here is the tunneling rate between the ith one- 

particle eigenlevel of the dot with spin a and the left 
(right) lead, Ei a is the one-particle eigenenergy of this 
level, P cq (N) denotes the equilibrium probability to find 
TV electrons in the dot (we assume the typical experimen- 
tal situation where the Coulomb blockade only allows N 
and N+l electrons in the quantum dot), P(Ei a \N) is the 
canonical probability to have the zth level of the spin-a 
sector occupied given the presence of N electrons in the 
dot, and f(E — fi) is the Fermi function at the effective 
chemical potential fj, which includes the charging energy. 



P eq (N)P(E irT \N)[l-f(E ia -^)}. (8) 



In Eq. 



is distributed according to the Porter- 



Thomas distribution for the GUE P 2 (T) = I exp(-r/r), 
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which only depends on the mean value T of the distribu- 
tion (we assume this mean value to be the same for the 
coupling to the left and right lead in the following). 

At zero temperature, only one level (ii,<7i) contributes 
in Eq. © such that fi = E ll(7l , P(E ia \N) = 1, and 
P eq (N) = 1/2. Thus, the zero temperature average con- 
ductance is given by (G) = j^-^fj- and the ratio of 
standard deviation to mean-value becomes a(G)/(G) = 
2/V5. Here, we have used (r£/(r£+r£)) = 1/3 and 
(r£ 2 /(Tf CT +r^) 2 ) = l/5 for the GUE distribution. At 
low temperatures, there are a few realizations of the 
RMT eigenlevel distribution where a second level (12,02) 
is within an interval of order k^T around the first level 
at the Fermi energy. Then, the second level also con- 
tributes to the conductance through the quantum dot. 
Neglecting the shift of the chemical potential due to the 
second level (i.e., keeping ^ = Ei 1(7l ), we calculated this 
two level situation. This gives the leading behavior in 
k B T/A for Eq. ©: 
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(9) 



For general temperatures, we maximize numerically 
the conductance Eq. (JSJ w.r.t. /i, and averaged over 
100000 RMT realizations of the eigenenergies and the 
dot-lead couplings. The results are shown in Fig. in 
comparison to the experiment of Patel et al. In contrast 
to the standard RMT result 10], the RMT ensemble for 
weak spin-orbit coupling describes the width of the con- 
ductance distribution and its change with temperature 
reasonably well at low temperatures, without any ad- 
justable parameter. Compared to the standard GUE, the 
width of the distribution is reduced at low temperatures 
because of the absence of level repulsion for levels with 
opposite spin. This results in higher probability to find a 
close-by level (with opposite spin and independent tun- 
neling rate), and leads to more RMT realizations in which 
two or more levels contribute to the low-temperature 
conductance. Having more independent channels for 
the conductance makes the probability distribution more 
Gaussian and decreases its width. At higher temper- 
atures, the experimental results are not adequately de- 
scribed b y sp in-orbit effects alone. In the regime k B T > A 
however [17]. one has to account for inelastic scatter- 
ing Tin. Taking the limit Ti n — > 00, we obtain the high 

temperature asymptotic behavior = J ^j^pp which 
gives reasonable results except for the quantum dot with 
diamond-symbols. Note that upon reducing fc^T/A, r; n 
will decrease, resulting in a crossover from the dot-dashed 
to the dashed line in Fig. [2] The inelastic scattering rates 
of (l?! would imply that the dashed Ti n = 00-line is ap- 
proached in the range fcsT=1.5-4 A. 

In Fig. we compare the full probability distribu- 
tion with the experimental one [10( at fcsT = 0.1A and 
fcgT = 0.5A. Within the experimental statistical fluc- 




FIG. 2: Width of the conductance distribution a(G)/{G) vs. 
temperature. At low temperatures, the RMT ensemble for 
weak spin-orbit interaction [dashed line; dotted line: low tem- 
perature behavior Eq. well describes the experiment |ic|| 
(symbols correspond to slightly differe nt q uantum dots), in 
contrast to standard RMT (solid line) 10]. At higher tem- 
peratures, a further suppression is due to inelastic scattering 
processes (dot-dashed line: Tin = 00 high-T asymptote). 



tuation, a good agreement is achieved without any free 
parameters, much better than for the standard RMT (lfij . 
This suggests that the spin-orbit strength is sufficient to 
fully decorrelated the spin-up and spin-down ensembles. 
With an estimate of the experimental Thouless conduc- 
tance (jT ~ 20 obtained from qt ~ V~N, this means that 
a spin-orbit coupling strength Qso ^ 10 -2 is required 
in the quantum dot of Ref. 10 (where we set \ = 1 in 
Eq. ©). In general, the crossover to the weak spin-orbit 
regime occurs at Qg (kpL) 5 ^ 2 > 1. Thus, the size of 
the dot and the parameter Qso which depends on the 
dot's specific asymmetry of the confining potential de- 
termine whether this quantum dot is in the weak spin- 
orbit regime. The size dependence might explain why 
earlier measurements by Chang et al. [Isj using very 
small quantum dots showed agreement with the standard 
RMT without spin-orbit interaction. A similar agree- 
ment was found by Folk et al. Il9 l , despite using similar 
large quantum dots as in Ref. llOt The contradictory re- 
sults of [19( and |l(| might be due to the better statistics 
of the latter experiment, or could be explained within the 
framework presented here, as following from differences 
in the confining potential (which might be, e.g., caused 
by differences in the realization of the two dimensional 
electron gas and the gate voltage), translating into dif- 
ferences in Qso- Alternatively, it is possible that the 
spin-orbit effect in both samples is weak, and the devia- 
tions from RMT in 10] should be explained by another 
mechanism (e.g. exchange [23^. 

In order to validate that the quantum dot is indeed 
in the weak spin-coupling regime described here, we sug- 
gest to repeat the experiment with an in-plane magnetic 
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FIG. 3: RMT predictions with weak spin-orbit coupling 
(dashed line) for the probability distribution of the Coulomb 
blockade peak conductance for a quantum dot at ksT = 0.1A 
(left figure) and fcsT — 0.5A (right figure), compared with 
the Patel et al. experiment [loll (histograms) and standard 
RMT theory (solid line). There are no free parameters in 
these distributions. 

field. A strong in-plane magnetic field should drive the 
system towards the strong spin-orbit scattering limit. In 
general, one would expect the spin-orbit scattering to 
suppress a{G)/(G). However, in the case of weak spin- 
orbit coupling the in-plane magnetic field, which drives 
the system towards a single GUE, regenerates the level 
repulsion. Therefore, we predict <j(G)/(G) to increase 
upon applying an in-plane magnetic field at low temper- 
atures. Another crucial test to the the weak-spin orbit 
scenario is the behavior in the absence of a magnetic 
field. In this case, the degeneracy is preserved but the 
spin-orbit coupling drives the system from the Gaussian 
orthogonal to the unitary ensemble Q. One implication 
is a strong suppression of the magnetoconductance. 

Finally we note that the d isag reements between RMT 
predictions and the results of |l0| can not be attributed to 
dephasing. Had this been the case, this experiment would 
indicate an appreciable dephasing even at low tempera- 
tures, in contradiction with theoretical predictions |2fj| . 
However, recent measurements of the low-temperatures 
dephasing rates [2l| are consistent with theory an d 
furthermore, it has been shown by Rupp and Alhassid 
22] that dephasing alone can not explain the results of 
10(. Our calculation shows that the spin-orbit coupling 
without dephasing can describe the low-temperature part 
of 01 > an d that the inclusion of strong dephasing gives 
reasonable agreement for the high-temperature part. 

In conclusion, we analyzed the effect of weak spin- 
orbit coupling on closed quantum dots in the presence 
of a perpendicular magnetic field which breaks the time- 
reversal symmetry. In this regime which can be real- 
ized for (some) quantum dots, the spin-orbit coupling 
does not lead to one non-degenerate GUE ensemble but 
to two independent GUEs for spin-up and -down elec- 
trons. This has important consequences, in particular, 
at low temperatures, as there is no level-repulsion for 
levels with opposite spins. The statistical distribution of 
the conductance peak maximum shows a good agreement 



with recent experimental distributions by Patel et al. llOj , 
but disagrees with experiments for similar sized quantum 
dots |l9j. The exchange interaction might yield similar 
changes in the statistical distribution |23j . and it is un- 
clear at present whether the complete explanation for 
the peak heights statistics behavior is given by the weak 
spin-orbit RMT. More experiments are needed to clarify 
the relative importance of the two effects and to explain 
the experimental contradiction mentioned above. If the 
spin-orbit effect is dominant, we predict an increase of the 
width of the distribution upon applying a strong in-plane 
magnetic field and a very low magnetoconductance. We 
further note that without spin-degeneracy there will also 
be no (5-function-like contribution in the level-spacing dis- 
tribution, in contrast to standard RMT. 
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